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ABSTRACT 

We investigate roles of the bulk comptonization process in the propagation of thermal photons 
emitted at the shock breakout of a supernova explosion. We use a hydrodynamical model based on a 
self-similar solution for the shock breakout. The propagation of photons is treated by a Monte-Carlo 
technique. Results of the simulations successfully reproduce the power-law spectrum observed by X- 
ray observations for XRF 080109/SN 2008D, if a mildly relativistic shock velocity is assumed. Effects 
of some radiative processes, photoionization, radiative recombination, and free-free absorption on the 
propagation of emitted photons are also investigated. It is found that none of them hardly changes 
the spectrum regarding the progenitor stars of type lb or Ic supernovae. Light curves calculated under 
the assumption of a spherical explosion indicate that the progenitor radius is required to be ^ 10^'^ 
cm. 

Subject headings: X-rays: bursts - shock waves - radiation mechanisms: non-thermal - supernovae: 
general - supernovae: individual (SN 2008D) 



1. INTRODUCTION 

The collapse of the iron core via photo-disintegration 
of nuclei at the final evolutionary phase of a massive 
star results in the formation of a strong shock wave in 
the deep interior. The shock wave propagates in the 
envelope, finally reaches the surface, and outshines the 
emission from the progenitor. This phenomenon occur- 
ring at the beginning of a supernova explosion is called 
supernova shock breakoutdKlein k Chevalieri ll978l: iFalkl 
119781 : iMatzner fc McKed Il999( ). A growing example 
of detections o f the shock break o uts, e.g., SN 2008D 
(iSoderberg et a l. 2008; ED l2008t iMazzali et all [20081: 
Modiaz et al. 200 9i), SN LS-04D2dc (iSchawinski et al ' 



2008; Gezaiiet_aL mM), SNLS-06Dlid (IGezari et a! 
20081 ). promotes investigations into its property, de- 
tectability, and usag e for cosmological studies (e.g. 
iTominaga etalMom . 

Among the example detections, the X-ray outburst 
XRF 080109 was serendipitously detected during the 
follow-up observations of a previously discov ered Type 
Ibc supernova, SN 2007uy (|Soderberg et al.ll200 8). It 
was associated with another Type Ib/c supernova 2008D. 
The parent galaxy was a spiral galaxy NGC 2770 at the 
distance of 28 Mpc. Analyses of the observed X-ray spec- 
trum and light curve have evoked a controversy on the 
origin of the emission. Some authors attribute the origin 
of the X-ray emission to the emergence of a radiation- 
dominated shock wave from the envelope of the progeni- 
tor, i.e., the supernova shock breakout. Others consider 
that SN 2008D is a gamma-ray burst (GRB) associated 
supernova and the X-ray emission is originated from the 
ej ected matter moving a t relativistic speeds. 

ISoderberg et al.l (|2008l l presented the X-ray light curve 
of XRF 080109 and results of the spectral analysis. The 
emission lasted for At ^ 600 s and its peak luminos- 
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ity was Lx = 6.1 x 10^"^ erg s^^. The total energy of 
the outburst is Ex = 2 x 10'*^ erg. The X-ray spectrum 
is best fitted by a single power-law with the photon in- 
dex ofF = 2.3±0.3. Assuming the presence of rela- 
tivistic ejecta, they estimated the bulk Lorentz factor at 
7 ~ 90. Because the scale length of the emission region, 
7^cAt ~ 10^'' cm, is much larger than the expected ra- 
dius of the progenitor star, R = 10^° cm, they excluded 
the possibility of the relativistic ejecta origi n of X-ra 
photons. Here c denotes the speed of light. Li (200^ 
performed a spectral analysis of XRF 080109 and found 
that the X-ray spectrum is well fitted by a power law 
spectrum with the photon index of F = 2.29^0 26 
sorbed by matter with the total equivalent hydrogen col- 
umn density of Nh = 6.83l} 3 x 10^^ cm~^, leading to a 
conclusion that the shock breakout origin is unlikely be- 
cause no black body component is found. IMazzali et al.l 
(|2008f ) also insisted on the relativistic ejecta origin be- 
cause of little contribution from a black body compo- 
nent to the X-ray spectrum. However, iModiaz et al.l 
(2009) claimed that the X-ray spectrum is well fitted 
by a black body combined with a power law spectrum 
with 50% bolometric contribution from the black body 
component. Their best fit values for the photon index 
F and the black body temperature Tbb are F 
and kTBB = 0.10 ± 0.01 keV, respectively. 

While X-ray observations have not led to agreement on 
the origin of the X-ray emission, observations in other 
wavelength ranges support the shock breakout origin. 
iBietenholz. Soderberg. fc Bartell ()2009f ) have performed 
very-long-baseline interferometry (VLBI) observations at 
30 and 133 days after the explosion. They succeeded in 
resolving the sizes of the source at both epochs and de- 
termined the 3cr upper limit to the expansion velocity 
to be ~ 0.75c , which rules out the presence of ultra- 
relativistic ejecta. From optical spectros copic observa- 
tions for SN 2008D in the nebular phase, iTanaka et al.l 
(|2009, ) found a double-peaked [O I] line profile, which is 
strong evidence of an aspherical explosion. They argued 
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that SN 2008D is a side- viewed bipolar explosion. 

If relativistic ejecta are responsible for the observed 
X-ray emission, X-ray photons will be emitted by syn- 
chrotron mechanism, which naturally explains the ob- 
served power-law (or black body + power law) X-ray 
spectrum. On the other hand, early investigations into 
the supernova shock breakout (iKlein fc Chevalie"il 119781 : 
iFalkl 119781: iMatzner fc McKed |1999D predicted thermal 
black body radiation for a UV/X-ray flash associated 
with the shock breakout. Therefore, some additional 
processes producing non-thermal photons are required to 
explain the observed power-law spectrum in the frame- 
work of the shock breakout origin. The bulk comptoniza- 
tion is a process to produce non-the rmal photons in the 
presen ce of a bulk flow of matter (jBlandford fc Pavni 
Il981air0 ). A small fraction of photons are repeatedly 
scattered by electrons across the shock front and grad- 
ually increasing their energies. The first application of 
th e process to t he su pernova shock breakout was done 
by iWang et"aLl ()2007D , who calculated the propagation 
of photons around the front of a radiation-dominated 
shock by a Monte-Carlo method. However, the shock 
profiles assumed in their calculations might oversimplify 
the situation and effects of other radiative processes com- 
peting the electron scattering, e.g., photoionization, ra- 
diative recombination, and free-free absorption, were not 
discussed. In addition, they did not construct a light 
curve from their model. Hence we try to e stablish a 
more realistic model than iWang et al.l (|2007| ) by using 
an appropriate hydrodynamical model, the self-similar 
solution for the supe rnova shock bre akout. The solution 
originally derived by iSakurail (|1960[) describes the prop- 
agation of a shock wave in a planar atmosphere with 
equations of non-relativistic fiuid dynamics. We treat 
the radiative transfer of thermal photons emitted by the 
shock front in the test-particle limit. Then we construct 
not only spectra but also X-ray light curves from mod- 
els with several different parameter sets. To derive the 
spectra, we neglect the contribution from photons inside 
the photosphere, which is defined as the position where 
the optical depth measured from the stellar surface be- 
comes unity. Using the results, we discuss the origin of 
the X-ray emission from XRF 080109. 

This paper is structured as follows. In §2, we describe 
formulation of the problem. The derivation of the self- 
similar solution for the shock breakout used in §2 is re- 
viewed in Appendix. An application of our model to 
XRF 080109/SN 2008D is presented in §3. Finally, §4 
concludes this paper. 

2. FORMULATION 

In this section, we explain qualitative natures of the 
supernova shock breakout, describe details of our model, 
and formulate the problem. 

The collapse of the iron core of a massive star results 
in the formation of a strong shock in the deep interior 
and the shock propagates toward the envelope. When 
propagating in the envelope, the shock is mediated by 
radiation. In the diffusion approximation, the average 
velocity Vdis of photons diffusing in a medium with the 
optical depth r is expressed by the speed of light divided 
by the optical depth, Vdis — c/t. Here r represents the 
optical depth from the outer edge of the envelope to the 
shock front. When the shock velocity 14 is larger than 



the diffusion velocity, photons behind the shock front are 
trapped . On the other hand, when the shock velocity 
is smaller than the diffusion velocity, photons can over- 
take the shock front and observers can see these photons 
as a precursor of the shock front. Therefore, when the 
relation Vs = Wdiff is satisfied, namely, the optical depth 
decreases to the ratio of the speed of light to the shock 
velocity, r — c/Vg, the shock breakout occurs. 

2.1. Shock profiles 

A semi-analytical description of time evolutions of hy- 
drodynamical variables at t he super n ova s hock break- 
out has been obtained by ISakurail ()1960[) . The so- 
lution was extended to the u ltrarelativistic case by 
iNakavama fc Shigevamal (|2005l ). Furthermore, the 
present authors extended the solution so that it can 
treat the energy deposition or loss at the shock front 
(jSuzuki fc Shigevamal 120071 ) ■ In this work we use the 
self-similar solution discovered by SakuraJ ()1960D to ob- 
tain velocity and density profiles behind the shock front. 
In the f o llowin g subsections, we review the approach of 
ISakurail ()1960t ) and describe the procedure to obtain the 
profiles. 

2.1.1. Self-similar solution 

ISakurail ()1960l ) studied the propagation of a strong 
shock in a medium whose density profile pq{x) can be 
written by a power of the distance x measured from the 
outer edge of the envelope of the progenitor star. 



kix^ 




for X > 0, 
for X < 



(1) 



where fci and S are constants characterizing the density 
structure of the envelope. Because we are concerned with 
the radiative envelope of a massive star, we take 5 = 3 as 
the fiducial value. The value of ki is treated later. Time 
evolutions of the velocity u{x, t), the density p{x, t), and 
the pressure p{x, t) behind the shock front are governed 
by the following equations for adiabatic gases in the plane 
parallel atmosphere. 
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where 7 represents the adiabatic index. In this paper, 7 is 
fixed to 4/3 because the shock is assumed to be radiation- 
dominated. The time t is measured from the point when 
the shock reaches the outer edge of the envelope, x = Q. 
For the self-similarity, we assume that the shock velocity 
Vs is expressed by a power of the position X of the shock. 



V. = koX 



-A 



(5) 



where k2 is a constant and we have introduced a param- 
eter A, the value of which is determined so as to satisfy 
appropriate boundary conditions. The procedure to de- 
termine A is summarized in Appendix. Integrating the 
definition of the velocity of the shock, dX/dt = 14, with 
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respect to time, the position of the shock is obtained as 
a function of time t. 



X = [fc2(l + A)(-t)]i/(i+^). 
Then we define the similarity variable ^ as 



A+l 



(6) 



(7) 



and introduce the non-dimensional velocity /(^), pres- 
sure g{£,), and density /i(^) as 



u{x,t)=k2X-^f{0, 

p{x,t) = hX'h{0. 



(8) 



Using these functions, the governing equations (HJ-Q 
are converted to ordinary differential equations. 



_,dln/i df S-X 



(9) 



A 



hd^ 1 + A 



9 S-2Xq , , 



Next, we consider the boundary condition at the shock 
front. The Rankine-Hugoniot relations for a strong shock 
determine the values of the velocity Uf, the pressure pf, 
and the density pf immediately behind the shock front, 



-Vs, Pi 



pi 



7 + 1 



7-1-1 °' 7+1'" " ' " 7-1 

The boundary condition at the shock front yields 



Po- (12) 



/(I) = 



1 



7 + 1' 



Ki) = 



7+1 



Mi) = 



7 + 1 
7^' 



(13) 



Integrating Eqs. ©-((Ill) from ^ = 1 to ^ = for the 
eigen value of A, we obtain the profiles of /(^), g(C)i ^'i^d 
/i(^). Otherwise the solution could not reach ^ = 0. 

2.1.2. Denvation of the shock profiles 

Once we obtain the profiles, /(C), ff(0, ^-^^d h{£), we 
can derive the time evolutions of hydrodynamical vari- 
ables u{x,t), p{x,t), and p{x,t) by using the relations 
© . To apply the solution to a specific problem, we need 
to fix the free parameters ki and ^2 introduced above. 

To fix them, we determine the position Xi and velocity 
Vi of the shock front at the shock breakout when the 
optical depth is equal to the ratio of the speed of light 
to the shock velocity. Since we can calculate the optical 
depth T by simply integrating the density profile ahead 
of the shock front, 



Xi 



■5+1 



kix dx = 



arkiX'. 
{S + l)pmH ' 



(14) 



where p{= 0.6) is the mean molecular weight, toh the 
Hydrogen mass, and tiT is the Thomson cross section, 
we obtain the value of ki for given Xi and Vi as. 



fci 



{5 + l)p,niHC 
<7TX'+'\Vi\ ■ 



(15) 
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Fig. 1. — Snapshots of the density profile (upper panel) and the 
velocity profile (lower panel). Each line represents the density or 
the velocity at t = —7.1 s (solid), —5.6 s (dashed), —4.1 s (dotted), 
and — 3.4 s (dash-dotted). 

The parameter k2 is obtained from the definition of the 
shock velocity ([S]), 



k2 = ViX^ 



(16) 



Next, we define a value X{ so that the optical depth 
T becomes unity when the shock reaches the position 
X = Xf. Using Eqs. (ITU) and (fT^ . the position X{ is 
expressed as 



l/(l+'5) 



Xi, 



(17) 



which corresponds to the final velocity V{ of the shock 
wave given by 



Vt 



m 

c 



-A/(l+5) 



Vi. 



(18) 



We stop the calculation for the photon propagation when 
X = Xf. 

An example of the thus derived profiles of the velocity 
u{x, t) and density p{x, t) are shown in Figure[TJ The free 
parameters are chosen as Xi ~ 10^^ cm and Vi = 0.3c. 

2.1.3. Uncertainty of the profiles 

Although the set of the derived profiles is an exact so- 
lution of the governing equations for non-relativistic adi- 
abatic gases, it partly deviates from a realistic model for 
a radiation-dominated shock propagating in the stellar 
envelope. Especially, the temperature profile after the 
shock breakout must deviate from that derived from the 
density and / or pressure profiles by making use of an ap- 
propriate equation of state because the present method 
does not take into account the energy transfer by photons 
overtaking the shock front. Therefore, in the following 
calculations of the photon propagation around the shock 
front, we only use two profiles, the velocity and density 
profiles. The former is used to determine the bulk ve- 
locity of electrons scattering photons and the latter is 
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used for the calculation of the optical depth. On the 
other hand, the pressure profile indicates a roughly con- 
stant post-shock pressure Ps — 10^°g/cni s^. Since the 
post-shock pressure is expected to be dominated by radi- 
ation pressure, Ps — aT^^^/3, where a is the radiation con- 
stant, we obtain the photon temperature of the order of 
fcTph ~ 0.1 keV, which is con sistent with results of spec- 
tral analysis of XRF 080109 (|Modiaz et al.l[2009l ). Then, 
we assume a uniform temperature profile with kTph = 0.1 
kcV behind the shock. 

2.2. Radiative processes 

Before we describe the procedure to deal with the pho- 
ton propagation, we estimate relative contributions of 
some radiative processes important in the stellar enve- 
lope. For photoionization and radiative recombination, 
we estimate the time scale of each process for oxygen, 
which is one of the most abundant heavy elements. For 
radiative recombination, we compare the absorption co- 
efficient with the scattering coefficient for Compton scat- 
tering. In the following discussion, we assume the enve- 
lope of a compact star such as Wolf-Rayet stars because 
the velocity of the shock wave propagating in the enve- 
lope is expected to reach mildly relativistic speeds, which 
is required for the bulk comptonization to work well. 

First of all, we estimate the average electron number 
density fie of the envelope when the electron scattering is 
a dominant source of opacity. Then, we suppose that the 
shock front is located at a distance AR from the outer 
edge of the envelope. Using the average number density 
fie, we can express the optical depth from the outer edge 
of the envelope to the shock front as nefx Ai?. Since the 
optical depth at the moment of the shock breakout is 
given by r — c/Vs, the average electron number density 
of the envelope is estimated as, 

arARVs ^10^ cm/ \0.'Sc J 

(19) 

Therefore, for a progenitor star with the radius R, the 
number of electrons ahead of the shock front is estimated 
as 

47rficAi?i?2=4 X 10"*^ ( ) 

V3xl0i5cm-V 

(lo^) (lo^) • ^^^^ 

2.2.1. Photoionization 

Next, we consider photoionization of heavy elements. 
X-ray photons emitted from the shock front ionize ele- 
ments in the envelope. Especially, heavy elements like 
oxygen and neon efficiently absorb the X-ray photons 
because the absorption edges of such elements are in the 
energy range of the X-ray photons. For example, the pho- 
toionization cross section of Hydrogen-like oxygen ion for 
photons with energy hv is given by 

(jRvbicki fc LightmanI 119851 ) . where 5(~ 1) is the gaunt 
factor, ais is the fine structure constant, and qq is the 



Bohr radius. The energy x is the ionization potential for 
oxygen, 

X = •yioF-m^c^ = 0.87 keV. (22) 

Therefore, for photons with the energy x, the photoion- 
ization cross section is of the order of dph — 10""^^ cm'^, 
which is much larger than the Thomson cross section ut ■ 
Namely, photoionization is a dominant opacity source 
until all the heavy elements become fully ionized. In 
fact, however, photoionization hardly deforms the X-ray 
spectrum because all the relevant elements become fully 
ionized immediately. To show this, we estimate the total 
emitted energy Ey^ of non-thermal X-ray photons. As- 
suming the photon temperature of Tph ~ 0.1 keV and the 
emitting radius of i? = 10^^ cm, the luminosity Lth of the 
thermal emission is obtained by the following formula, 

L.,. . 4rf.se3?„ ^ 10« (^)^ (^)*erg/s. 

(23) 

where (tsb is the Stefan-Boltzmann constant. Since the 
time required for the shock wave to propagate across a 
shell with the width Ai? is given by AiJ/Fg, the total 
energy £^th of the thermal emission is 

Ai? 42 R Y f Tpi-, Y 

Furthermore, introducing a parameter e, we express the 
total emitted energy of non-thermal X-ray photons as 
Ex — eEth ~ lO'^^e erg. Once the outburst energy Ex 
is given, we can estimate the number of X-ray photons 
having the energy of a few keV and the number is of the 
order of lO^^e photons, which is larger than the number of 
electrons estimated above (l20l) when a substantial con- 
tribution of non-thermal photons, e > 0.1, is assumed. 
From the estimated value of the number of photons and 
the volume of the envelope, the number density yields 

lO^^e 
"P^"47ri?2Ai? 

= 8xl0-m f^y\m-,(25) 
\0.lJ Vo.lkeVy \Q.3cJ ^ 

Using these values, we can evaluate the time scale for 
photoionization as 

Tph = i = 4xlO-3s, (26) 

ctTphriph 

which is much shorter than the time scale for the shock 
wave to propagate in the envelope. 

2.2.2. Radiative recombination 

Then, we estimate the time scale of radiative recombi- 
nation for fully ionized oxygen. The time scale of radia- 
tive recombination for ions with the atomic number Z is 
expressed by 
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The radiative recombination coefficient a^'^ is given by 

X (0.4288 + 0.5 In /3 + 0.469/3"^/^) cm^s, (28) 
(|Seatonlll959|) . Here the parameter (3 is given by, 

(29) 



^ fcT V O.lkeV 



Using the above relations and the temperature of kT = 
O.f keV inferred from the spectral analysis of the X-ray 
spectrum (Modjaz e t al. 2009), we obtain 



Trad = 3 X lO-^i;- 



oxygcn 1 3 ^ -^Qi5 



S, 



(30) 



where l^xygcn is the number of oxygen atoms per nucleon. 
Since the estimated time scale for radiative recombina- 
tion is much longer than that of photoionization even for 
pure-0 atmosphere (^oxygen = 1), we can safely assume 
f ully ioniz e d mat ter in the envelope. 

iBadnelll ()2006D provides results of the calculation of 
dielectric recombination rate coefficients for hydrogen- 
like ions. His results show that the dielectric recombina- 
tion rate coefficient for the fully ionized oxygen is smaller 
than the radiative one for the temperature range of our 
interest. Therefore, we can neglect the contribution of 
dielectric recombination. 

2.2.3. Free-free absorption 

Finally, we consider free-free absorption. The free-free 
absorption coefficient for ions with the charge Z is given 
by 

= 3.7 X 108T-i/2z2nenii.-3(l - e-''''/'=^)gff cm-\ 

(31) 

(e.g. iRvbicki fc Lightmanl ll985t). where rn is the num- 
ber densities of ions, ly is the frequency of the incident 
photon, and gg represents the gaunt factor. In the fol- 
lowing, we consider the contribution from oxygens. For 
the temperature deriv ed by the spectral analysis of XRF 
080109, kT = 0.1 keV (Mazzah et al."2008': 'M odjaz et al.l 
[2009^ . the free- free absorption coefficient is estimated to 
be aJJ ~ 10~^^cm~^ for hiy — 0.1 keV, which is much 
smaller than the scattering coefficient for Compton scat- 
tering, fieCTT — 10"'^"'^ cm""'^. Here we have assumed that 
the gas is fully ionized and that the gaunt factor is unity. 
On the other hand, for the surface temperature of a mas- 
sive star, T ~ 10** K, the free-free absorption coefhcicnt 
for low enerra^ photons is comparable to the scattering 
coefhcicnt, ~ 10~^°cm^^ for hu = 10 eV. However, 
since the effective optical thickness r* of the envelope, 
which is given by 



(32) 



(e.g. IRvbicki fc Lightmanl 119851 ) is smaller than unity, 
Th. ~ 0.3, even for T ~ lO'* K and hv = 10 eV, most 
photons must escape from the envelope without being 
absorbed. As discussed in ^12.2.11 matter in front of the 
shock front is exposed by the X-ray photons and free elec- 
trons thereof must have higher temperatures. Hence we 
may overestimate the effective optical thickness. For the 



progenitor stars of type lb or Ic SNe, the contribution 
from other elements, e.g., helium and carbon, may be 
dominant. However, the difference of the composition of 
the envelope, does not alter the conclusion, r* < 1. Thus, 
we neglect effects of free-free absorption in the following 
calculations. 

2.3. Radiative transfer 

In this subsection, we describe our method for com- 
puting the generation and the propagation of photons in 
given velocity and density profiles. 

2.3.1. Photon generation 

First of all, we generate seed photons having a black 
body spectrum with the photon temperature Tph as- 
sumed to be fcTph = O.lkeV. The frequency v of each 
photon is determined by using random numbers so that 
the distribution f^{v) of photons obeys the following 
form, 

,,2 



exp(/iiy/fcrph) 



1 



(33) 



The direction vector of each photon is specified by 
the inclination angle 6 and the azimuth angle (p as 
(sin6' cos sin0 sin cos 0). For seed photons, the an- 
gles and (j) are determined so that the angular distri- 
bution of photons becomes isotropic. In our model, the 
seed photons are assumed to be produced at the shock 
front at a constant rate. Numerically, a fixed number 
of photons are injected at the shock front at every time 
step. 

For the velocity and density profiles around the shock 
front, we divide the spatial interval [0,1/] into equal 
zones with the width Ax = LjN^- Thus the spatial 
coordinate of the center of the z-th zone is evaluated as 



Ax(i - 1/2) for \ <i<N^ 



(34) 



We can derive these profiles at any time from the self- 
similar solution in §2.1 evaluated at each of xi. The 
interval must cover the region from the initial to the final 
position of the shock. In the following calculation, we use 
the mesh covering the space [0, 2X\\ (L — 2X\) with 1000 
zones [Nx = 1000). The total number of seed photons is 
about 10^ 

2.3.2. Photon propagation 

We calculate the path through which a photon prop- 
agates during a time interval At ~ — using an 
optical depth t" defined as 



(35) 



Here we approximate the region around the photon as a 
uniform shell with the electron number density n" . This 
density is evaluated at the position x where the photon 
is located at the time t — as 



{x - Xi)n2^,+i + [xi+i - x)n^i 



(36) 



Here the photon is supposed to be located in a position 
between Xi and Xi+i. 

The photon may be scattered by an electron within 
this time interval. The probability P of this event is 
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expressed by the optical depth as 

P = l-exp(-r"). 



(37) 



Generating a random number Ri ranging from zero to 
unity, it is determined whether the photon travels freely 
or is scattered by an electron. When the random number 
i?i is larger than the probability P, the photon travels 
freely. On the other hand, when Ri is smaller than P, 
then the photon travels by a distance of 



I 



ln(l - i?i) 



(38) 



before the scattering. In this case, we carry out the 
procedure for the scattering given in §2.3.31 After that, 
we evaluate the optical depth from the residual distance 
cAt — I as T = aTricicAt — I) and then repeat the pro- 
cedure above. 

2.3.3. Compton scattering 

Since there is a well developed Monte-Carlo tech- 
nique to d eal with Compton scattering of photons (e.g. 
[Pozdnva kov et al.| [l977. 1983,), we adopted this method 
and review the outline in this section. 

At first, we specify the momentum (px,Py,Pz) of an 
electron. In our model, the momentum distribution 
/c {Px TPy,Pz) of electrons at the scattering site is assumed 
to be a shifted Maxwell-Boltzmann distribution with the 
temperature To, 



fc{Px-,Py,Pz) oc exp 



{Px 



■Pz 



2mr.kTr. 



, (39) 



where me is the electron mass. In the following, the 
electron temperature Te is assumed to be the same as 
the photon temperature, i.e., kT^ — kTph — 0.1 keV. The 
other parameter characterizing the above distribution is 
the bulk velocity Uc of electrons at the scattering site. 
In the same way as we calculate the electron number 
density in i j2.3.2[ a linear interpolation of the velocity 
profile is used to calculate the value of Uc- For given 
Tc and Uc, the momentum of the electron is determined 
using random numbers. 

We express the 4-momentum of the photon before and 
after scattering in the rest frame of the electron as 



Pi = 



Ci/C 

e-i sin 9i cos ( 
Ei sin 6i sin (, 
ej cos 9i 



Pf = 



ef/c 
et sin 6{ cos < 
ef sin 6{ sin < 
ef cos 9{ 



(40) 

respectively. We use the Thomson cross section rather 
than the Klein-Nishina cross section for Compton scat- 
tering, because energies of photons considered here are 
significantly lower than the electron rest energy. Thus 
the differential cross section da/dVl is given by 



-crT(l 



cos 



0), 



(41) 



dfj 

(see e.g.. iRvbicki fc Lightmanlll985f) . Here we introduce 
the angle 8 between the directions of a photon before 
and after scattering, 

cos 9 = cos 9i cos 9i + sin 9i sin 9i cos((/)i — 4>f). (42) 



observer 




photoshere 



atmosphere ' 

Fig. 2. — A schematic view of the geometry of the emission region 
of X-ray photons. 

We determine the angles 9i and after scattering by 
two random numbers so that the probability of the pho- 
ton scattered into each solid angle follows the above dif- 
ferential cross section (|4ip . We obtain the energy of 
the scattered photon from the well-known formula (e.g., 
IRvbicki fc LightmanlfT985[) . 



TOcC 



mcC^ + ei(l — cos 8) 



(43) 



After we obtain the 4-momentum of the photon after 
scattering, we apply the inverse Lorentz transformation 
to the 4-momentum pf to obtain that in the original 
frame. 

2.4. Derivation of light curves 

Since the adopted self-similar solution assumes the 
plane-parallel atmosphere and we only consider a region 
near the surface, we need to take into account several 
effects when deriving light curves of photons produced 
by the shock wave. We describe the procedure to derive 
light curves in this section. 

After we stop the calculation of the propagation of 
photons around the shock front, we evaluate the opti- 
cal depth from the location of the photon to the surface 
Tph. First of all, we assume that non-thermal photons in- 
side the photosphere are thermalized before they escape 
from the stellar surface. Thus, we only deal with pho- 
tons whose optical depth is smaller than unity, Tph < 1, 
to derive light curves of non-thermal emission. Next, we 
calculate the propagation of the thus selected photons 
toward the surface according to the procedure described 
in §2.3. Then, for each photon, we obtain the angle 9ph 
between the direction of the propagation of the photon 
and the shock normal when the photon reaches the sur- 
face. Finally, we consider light travel effects by using 
the angle 0ph. The situation is schematically shown in 
Figure [21 We assume a spherical explosion. As a conse- 
quence, the shock breakout occurs simultaneously and in 
the same way at every points on the surface of the pro- 
genitor star. Then, we assume that a photon with 9ph is 
emitted at a point on the surface of the star, where the 
angle between the line of sight and the radial direction 
is equal to ^ph- The time t is measured from the point 
when photons with 0ph = are observed. Then photons 
with 9ph are observed at t = R{1 — cos9ph)/c due to a 
difference of the path from the emitting point (see Fig- 
ure [2]), where R is the radius of the photosphere of the 
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photon energy [eV] 

Fig. 3. — An example of the resultant energy spectra of photons 
experienced the bulk comptonization process in supernova shock 
breakout (solid line). The free parameters characterizing the shock 
profile are Vi = 0.3c and Xi = lO^^cm. The radiative envelope of 
a massive star (5 = 3) is assumed. A black body spectrum with 
a photon temperature of 0.1 keV (dashed line) and a power-law 
spectrum with the exponent of —2.4 (dotted line) are also plotted. 

star. Taking the above two effects, i.e., the optical depth 
and the light travel effect, into account, we derive light 
curves of X-ray photons emitted after supernova shock 
breakout. 

3. RESULTS 

We calculate spectra and light curves of photons hav- 
ing experienced the bulk comptonization after supernova 
shock breakout according to the procedure described in 
the previous section. In this section, we apply our model 
to the observed X-ray spectrum and light curve of XRF 
080109/SN 2008D. 

The progenitors of SNe lb and Ic are expected to be 
Hydrogen-deficit compact stars, whose typical radius is 
of the order of 10^^ cm. On the other hand, the long 
duration of XRF 080109, At ~ 600 s, predicts a much 
larger scale of the emitting regio n, i.e., cAt ~ 2 x 10^^ 
cm. However, as previous works (ICampana et al.l 120061 : 
iWaxman et al.l 120071: iWang et al.l I2007D have claimed, 
we can attribute the contradiction to the presence of a 
dense stellar wind. Namely, the shock breakout occurs 
in the dense wind that has the photosphere with a much 
larger radius than the stellar envelope. Then, we adopt 
R = 10^^ cm as the radius of the photosphere. This 
modification does not alter the conclusion in i.e., 
we can neglect photoionization, radiative recombination, 
and free-free absorption. 

3.1. Spectra 

Figure [3] shows an example of energy spectra of pho- 
tons calculated from a shock breakout model presented in 
Figure [TJ To obtain the spectra, we use photons outside 
the photosphere (r = 1). The values of the free parame- 
ters taken for this model are Vi = 0.3c and = 10^^ cm. 
The dashed line represents a black body spectrum with 
the photon temperature of fcTph = 0.1 keV. While the re- 
sultant spectrum is consistent with the black body spec- 
trum in a low energy range (10 eV-1 keV), it significantly 
deviates from the black body above 1 keV. The dotted 
line represents a power-law spectrum with the exponent 



of —2.4, whic h is implied by the observ ed spectrum of 
XRF 080109 (ISoderberg et al.l [20081: lUI 12008'). The re- 
sultant spectrum above 1 keV remarkably reproduces the 
observed power-law. 

In order to check the validity of the calculation, we 
have calculated the energy spectrum of photons injected 
at a point which is in the downstream and is sufficiently 
far from the shock front. We confirmed that the spec- 
trum remains black body one with Tph = 0.1 keV. There- 
fore, we can attribute the formation of the power-law tail 
in the energy spectrum to the interaction of photons with 
the shock front. 

Figure |3] shows energy spectra of photons calculated 
for various Vi and fixed Xi{= 10^^ cm). Four different 
lines represent models with Vi — 0.1c (solid), V = 0.2c 
(dashed), V = 0.3c (dotted), and V = 0.4c (dash- 
dotted). We can easily recognize a tendency in the spec- 
tra that a larger shock velocity model has a more promi- 
nent high energy tail deviating from the black body spec- 
trum. While the power-law tail is absent in the spectrum 
of models with V < 0.1c, the exponent of the power- 
law part of the spectra converges to F ~ — 2 for models 
with the initial shock velocity larger than 0.4c. However, 
since the shock profile we used here are valid for non- 
relativistic shock velocity, the value of the photon index 
may change for relativistic cases. 

Finally, we investigate the dependence on the density 
profile. The density profile of the stellar envelope is char- 
acterized by the parameter S in our model. We performed 
calculation of models with 6 = 2 and S = A. The resul- 
tant energy spectra are shown in Figure[S] We found that 
these energy spectra also show power-law tails with the 
same index as the case for (5 = 3, which is a physically 
natural consequence. High energy photons comprising 
the power-law t ail are produced by repeat ed crossing of 
the shock front (jBlandford fc PavneliTOSTal lbl) . The prob- 
ability that photons diffusing out from the shock front 
scattered back to the downstream is determined by the 
optical depth of the envelope, which is given by integrat- 
ing the density profile, rather than the density profile 
itself. In addition, as long as the optical depth is larger 
than unity, most of photons must be scattered back to the 
downstream. Slight differences of the density structure 
do not alter the bulk comptonization process. Therefore, 
the results weakly depend on the parameter d. 

3.2. Light curves 

Using the procedure described in §2.4| we calculate 
light curves of X-ray photons produced at the shock 
front. One of the results is shown in Figure [SI The 
radius of the progenitor star is assumed to be 10^'^ cm, 
which is required to explain the observed duration of the 
outburst of a few hundred seconds. 

We mention an inadequacy of our light curve model. 
Results of optical spectros copic observat i ons o f the neb- 
ular phase of SN 2008D (jTanaka et al.l I20090 strongly 
indicate the aspherical nature of SN 2008D. On the con- 
trary, our model assumes a spherical explosion, which 
leads to a consequence that the shock breakout simul- 
taneously occurs at every point on the surface of the 
progenitor star. For an aspherical explosion, the shock 
breakout does not occur simultaneously in different ra- 
dial directions. The assumption of spherical symmetry 
is not appropriate when the duration of the shock break- 
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Fig. 4. — The energy spectra for various initial velocity Vi. Each 
line represents the model with V; = 0.1c (solid), V"; = 0.2c (dashed), 
Vi = 0.3c (dotted), and = 0.4c (dash-dotted). The initial posi- 
tion of the shock is fixed (Jf; = 10^^ cm). The radiative envelope 
of a massive star (S = 3) is assumed. 
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Fig. 6. — An example of light curves of X-ray photons (1 
keV< hi/ < 10 keV) calculated using the same model as in Fig. 
\3\ Different lines represent the numbers of photons in the energy 
ranges of 1-10 keV (solid), 1-3 keV (dashed), and 3-10 keV (dotted) 
as functions of time. 




TO 100 1000 lO" 

photon energy [eV] 



Fig. 5. — The energy spectra for various 5. Each line represents 
the model with <5 = 2 (solid) and 5 = 4 (dashed). The initial 
velocity and the initial position of the shock are fixed (Vi = 0.3c 
and Xi = 10" cm). 



out is comparable to that of the X-ray outburst. This 
effect modifies results of our light curve model. How- 
ever, the investigation into this asynchronism requires 
two-dimensional hydrodynamical simulations of super- 
nova explosions. Therefore, we consider that the model- 
ing of light curves of X-ray photons from an aspherical 
explosion is beyond the scope of this work. Furthermore, 
the light curve of X-ray emission from XRF 080109 have 
large errors due to the paucity of photons. Therefore, it 
is difficult to judge whether SN 2008D was a spherical 
or aspherical explosion using only our results assuming 
plane-parallel atmosphere. 

3.3. Energetics 

Finally, we check the consistency of our model with the 
observational results. Since we have assumed the black 
body radiation with the photon temperature of Tph = 
0.1 keV and the scale of the emitting region to be i? = 
10^'^ cm, the luminosity Lth of the emission is of the 



order of 10**^ erg/s. The time during which the shock 
wave emits thermal photons, {Xi — X{) /Vs is of the order 
of seconds. Therefore, the total energy of the thermal 
emission is evaluated as Eth = 10'^^ erg. On the other 
hand, from the results of the simulation (Fig. [3]), the 
ratio of the energy of non-thermal X-ray photons with 
0.3-10 keV to that of the thermal emission is evaluated 
as e = 0.3-0.4. Using these values, we obtain the total 
emitted energy of non-thermal X-ray photons with 0.3- 
10 keV, Ex = eEth =(3-4)xl0''s erg, which agrees with 
the observed value very well. In other words, adopting 
the progenitor radius of lO^'^ cm, we can explain both 
the duration and the total energy of the observed X-ray 
emission. It should be noted that the above estimate 
has an uncertainty. To obtain the spectrum shown in 
Fig. |3l we use photons outside the photosphere (r = 1). 
However, photons in the downstream of the shock must 
diffuse out later and contribute to the thermal part of 
the spectrum. Especially, photons in the region where 
T < c/Vs can escape from the shock front. Therefore, we 
may overestimate the ratio e by a factor of c/Vg. 

4. CONCLUSION 

In this paper, we calculate the propagation of photons 
at the supernova shock breakout by using a Monte-Carlo 
technique in order to explain the power-law spectrum re- 
vealed by X-ray observations for XRF 080109/SN 2008D 
in the framework of the shock breakout origin of the X- 
ray photons. 

We use the self-similar solution for the propa gation of 
a shock wave in a stellar envelope presented in iSakurail 
fil96Ql to acquire the velocity profile and evaluate the 
optical depth for photons around the shock front. Our 
hy drodynamical mod el is more realistic than that used 
in iWang et al.l ()2007[ ) and allows us to obtain plausible 
constraints on the nature of the shock breakout occurred 
in SN 2008D. Our simulations clearly demonstrate that 
a mildly relativistic shock velocity, Vg > 0.3c, is required 
for the spectrum of photons emitted from the shock front 
to deviate from the black body spectrum due to repeated 
scatterings around the shock, i.e., the bulk comptoniza- 
tion. The resultant spectra successfully reproduce the 
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power- law spectrum with the exponent of ~ 2.4 seen in 
the X-ray spectrum of XRF 080109. We also investi- 
gate whether some radiative processes, photoionization, 
radiative recombination, and free-free absorption, inter- 
rupt the formation of the power-law tail and find that 
none of them hardly changes the shape of the spectra as 
long as progenitor stars of type lb or Ic supernovae are 
concerned. 

Assuming a spherical explosion, light curves of the X- 
ray photons are derived. To account for the duration 
of the X-ray emission from XRF 080109, the progenitor 
radius of ~ lO^^cm is needed. Although this value is 
somewhat large compared to the typical radius of Wolf- 
Rayet stars, the shock breakout taking place in a dense 
stellar wind may explain the discrepancy. Furthermore, 
our model successfully explain the total energy of the ob- 
served X-ray emission. Because the density structure of 
the wind is not known, it is unclear whether our power- 
law approximation of the atmosphere is appropriate for 
the situation. However, the fact that our results weakly 
depend on the parameter 6 suggests that the bulk comp- 
tonization can operate even in the dense wind and pro- 
duce a power-law X-ray spectrum. 

Our model does not include the contribution of pho- 
tons from the downstream of the shock and the matter 
ejected at t > 0. However, as we noted in §3.1, photons 
originating from far downstream have a black body spec- 



trum. So they do not contribute to non-thermal photons. 
In the same way, photons emitted by the ejected matter 
are not important. In particular, after the shock break- 
out, the ejected matter cools rapidly due to adiabatic ex- 
pansion. It must emit optical or UV photons rather than 
X-ray photons. Therefore, both of them hardly change 
at least the high energy part of the energy spectra cal- 
culated using our model. As we noted above, however, 
these thermal photons may decrease the ratio of thermal 
to non-thermal emission from the value obtained in §3.3. 

It should be noted that our model to calculate the radi- 
ation from the supernova shock breakout may have a few 
inadequacies. The used self-similar solution cannot deal 
with highly relativistic flows and the photon temperature 
is assumed to be equal to the gas temperature even after 
the shock breakout occurs. Such caveats are improved 
by calculating the radiation from the supernova shock 
breakout by using more sophisticated simulation codes, 
e.g, a multi-group radiation hydrodynamics code. 
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APPENDIX 

DERIVATION OF THE SELF-SIMILAR SOLUTION 
In this section, we review the procedure to determine the eigen value A. The following steps are originally presented 



in Sakurai (196y). 
At first, we solve Eq. ([TT]) with respect to dg/d^ 



dg 
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df <5-(7 + 2)A - 
'^d^ 1 + A 



d^ 1 - ^/ 

Substituting the above expression into Eq. (|10p , the following equation is obtained, 

^9 



dC HI - C/) 



■ df ^-(7 + 2) A - 



A ^^_5-2Xg 



1 + A 



1 + A /i 



Next, subtraction of both sides of Eq. © from those of Eq. ([TU|) leads to the following equation. 



(1-e/) 



ding dlnh 



di di 

Here we introduce new variables y and z as functions of i^. 
Using these functions, Eqs. (jA2l) and (IA3p are transformed into. 



7^^ z2 ' 



(1-y) 



dlnz 
d\n^ 



and 



dlny 



(5- (7 + 2)A 
(1 + A)7 



y(z-l) + 



1 + A 



{z-rf 



a-2X 
(1 + A)7 



y, 



dlnz 2A~ (7- 1) 7 + 1 

d\n^~^^^ 'd\n^^ 1 + A ^+1 + A- 

Then, the elimination of the variable ^ from the above two expressions yields 

dz 7(z - l)(Az + 1) + [(2A - (5 - 7)z + 7] y 



dy 27 - 7(7 - 1) A - [,5 + (<5 + 2)7 - 2A]y + 7(7 + 1) Az ' 



(Al) 

(A2) 

(A3) 

(A4) 

(A5) 
(A6) 

(A7) 
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We solve this equation numerically to obtain the eigen value A. 

Here we consider the boundary conditions of the equation. From the definition of the variables (|A4[) and the boundary 
condition at the shock front (|13p . their values yt and Z{ at the front are obtained as 

Vi = 7, zf = A8 

7 — 1 7 — 1 

On the other hand, another boundary condition is obtained as follows. In fact, Eq. (|A7|) have a singular point. From 
Eq. (jA5p . we can easily find that the singular point is located at j/s — 1, where the L.H.S of Eq. (|A5[) vanishes. Then, 
a requirement that the R.H.S of Eq. (|A5p also vanishes at j/s = 1 leads to the value of z at the singular point, 

7A 



From the thus derived boundary values, we can determine the eigen value A so that Eq. (jA7l) can be integrated from 
{yi,Z{) to (j/sj-Zs). For 6 = 3 and 7 = 4/3, which are used in this work, the eigen value is found to be A = 0.5 5727. 
The values of A for various S and 7 are found in the original work (Sakurai 1960) or lSuzuki &: Shigevamal (|2007l) . We 
obtain the profiles /, g, and h by numerically integrating Eqs. (l^- (|lip for the eigen value. 
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